Exponents and bounds for uniform spanning trees in d dimensions 
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Uniform spanning trees are a statistical model obtained by taking the set of all spanning trees on 
a given graph (such as a portion of a cubic lattice in d dimensions), with equal probability for each 
distinct tree. Some properties of such trees can be obtained in terms of the Laplacian matrix on the 
graph, by using Grassmann integrals. We use this to obtain exact exponents that bound those for 
the power-law decay of the probability that k distinct branches of the tree pass close to each of two 
distinct points, as the size of the lattice tends to infinity. 



In graph theory, a tree on a graph is a connected subset 
of the vertices and edges without cycles, and a spanning 
tree is a tree that includes all n vertices of the graph (it 
must then have n — \ edges). Results for the number 
of spanning trees on a given graph go back to the nine- 
teenth century (see e.g. Ref. [1]). If each spanning tree 
is given equal probability, we obtain uniform spanning 
trees. In this paper, we consider uniform spanning trees 
on (a portion of) the square, cubic, or hypercubic lattice 
in d dimensions. One would like to characterize the frac- 
tal properties of the trees as the size (number of vertices) 
of the lattice goes to infinity. One characteristic is the 
probability that two well-separated points are nearly con- 
nected by fc = 1, 2, 3, . . . distinct branches of the tree, or 
alternatively by distinct paths along the tree, and these 
are expected to behave as power laws that are described 
by critical exponents. We will study these by methods 
based on the classical results, and obtain some exact ex- 
ponents, which serve as bounds for more general ones. 
(In two dimensions, the exact results have been known 
for some time [2-6].) The motivation to consider this 
problem came from its connection to some optimization 
problems, which are in turn connected with the ground 
states of classical systems with quenched disorder, such 
as Ising spin glasses. In the two-dimensional case, there 
is also a connection with loop models, the Q ^ Potts 
model, and Coulomb gases in conformal field theory [2-4] . 

First we note that the result [variously attributed ei- 
ther to Kirchhoff (1847), or to Sylvester (1857), Bor- 
chardt (1860), and Cayley (1856)] for the number TV of 
spanning trees on a graph can be written in the following 
generalized form: 



7V = cof A("^1'^i) = (_l)^i+!/i det A^^^i'^^i), 



(1) 



where we have recalled the definition of the cofactor. 
Here Xi, yi = 1, 2, ...label the vertices in the graph, 
the matrix A (the lattice Laplacian) is defined as 

{degx if X — y, 
—t if X and y are connected by t edges, 
otherwise, 

(2) 

and A^^^'^i^ means the minor of {xi,yi), i.e. A with row 
xi and column j/i deleted. For xi =2/1, this reduces to 



the better known result. The effect of deleting a row and 
column is to remove the zero mode that would otherwise 
cause the determinant of A to vanish. 

The result generalizes further to a relation that in- 
volves the number J\f(^iv^'^2V2^---'^kyk) spanning sub- 
graphs without circuits with k components, and with Xi, 
yi in the same component for each i (we will assume that 
all Xi, yi are distinct). The result is 



cof A(^i 
± 



■Xk-Vl-'-Vk) — 



E 



Here the cofactor is again 

(_l)i:ti(-.+y.)detA( 



xi...xk,yi---yk) 



(3) 



(4) 



where rows Xi and columns yi have been deleted, and P 
runs over permutations of k symbols. The overall sign 
on the right hand side depends on the details of how 
the vertices are labelled and is uninteresting. Both the 
generalizations are mentioned by Ivashkevich [5] (see also 
Ref. [7]), but he omits the signs in the cofactors. The 
results can be proved by an extension of the proof given 
for example in Ref. [1]. 

In the following we will consider a graph that is a 
bounded portion A of the d-dimensional cubic lattice Z"* 
(with edges that connect only nearest neighbors at Eu- 
clidean distance 1 in lattice units). We will be interested 
in the following property of a spanning tree. We choose 
two vertices x, y, together with a neighborhood of each. 
We assume that the neighborhoods are chosen in such 
a way that the boundary passes through some vertices, 
but no edges of A cross the boundary; all edges are either 
inside or outside. We take k vertices Xi on the boundary 
of the neighborhood x, and k vertices yi on the boundary 
of that of y. In practise, this can be satisfied using neigh- 
borhoods that are approximately balls of radius of order 
k^/id-i)^ We can now look at the part of the tree lying 
outside the two neighborhoods; this amounts to a forest 
of trees, with each tree rooted on both the boundaries 
of the neighborhoods of x and y. We ask whether, for 
each i, the points Xi, yi lie in the same connected com- 
ponent in this forest, and are in a distinct component 
from any other pair Xj, yj. If so, then in terms of the 
original tree the XiS are connected to the corresponding 
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FIG. 1. Example of a spanning tree on a portion of the 
square lattice, with points x, y and a neighborhood of each 
marked, with vertices xi, X2, xs and j/i, 2/2, ys on the bound- 
ary of each shown. In this example, there are paths, lying 
outside the neighborhoods, that connect Xi to for each i, as 
required, and those paths are distinct, but there arc only two 
distinct branches (connected components) outside the neigh- 
borhoods. The pairs X2, 2/2 and xs, ys do not lie on distinct 
branches. 

yiS by branches of the tree that are distinct outside the 

two neighborhoods selected, and we call this "crossing k 
times (between specified points in the neighborhoods of 
X and y) by distinct branches". See Fig. 1. For given A, 

neighborhoods of x and y, and points Xi, yi, denote the 



number of such trees AAi 



ixiyi,X2y2,---,XkVk) 



branches 



Then we also 



define corresponding probabilities, 

q~,(xiyi,X2y2,---,XkVk) _ \r{x\Vi,X2V2,---,Xkyk) i \r /r\ 
' branches ~ "'^branches I ■'^ ' 

It is obvious that there is a relation between this def- 
inition, and the fc-component spanning subgraphs with- 
out circuits considered before, if the latter are defined on 
A~, that is A with the interiors of the neighborhoods of 
X and y removed. Given a fc-component spanning sub- 
graph of A~, a spanning tree of A can be obtained by 
adding k — 1 edges, each inside either the neighborhood 
of X or of y, together with one edge for each vertex in 
A — A~ . If we also assume that there are exactly k ver- 
tices on the boundaries of each of the two neighborhoods, 
then is obtained by summing over all 

the ways to add edges to each fc-componcnt spanning 
subgraph of A~ (with Xj and y^ in the ith component 
for all i) to obtain a spanning tree of A whose branches 
are distinct outside the two neighborhoods, and then 
also summing over the /c-component subgraphs used in 
this construction. As the maximum possible number of 
ways to add edges is limited by the size of the neigh- 
borhoods used, AA^^anches^'''''''^'^''' °^^y Sightly larger 
than the number Ar(^iyi.a:2y2,...,xfei/fe) fc-component 

spanning graphs of A^. (More precisely, there is a bound 

^ranchef" ' '"""'^ < Cfe,d^(-i^i'-^^^. -'-'=^'-) where Cfc.rf is 

independent of the distance from x to y.) 



Another possible definition of crossing from x to y 
would require only that the k crossings of the spanning 
tree be made by k paths on the tree that are distinct 
(have no edges in common) outside the neighborhoods 
(for k = 1 these are the same thing, and the probability is 

one). See Fig. 1. These numbers A^^X'"''"""''"''"^ can 
be obtained by adding edges to the fc-componcnt span- 
ning subgraph of A~ in arbitrary positions in A, and so 

j^j-(xiyi,x2y2,---,XkVk) ^ j^{xiyi,X2y2,--;Xkyk) g^^^ 



paths 

Krixiyi,X2y2,...,Xkyk) I Kf < n:){xiyi.x2y2..--,Xkyk) 
•^^ — ' branches 

^ ^{xiyi,x2y2,---,Xkyk) 



(6) 



These contimie to hold even if there arc more than k 
vertices on the boundary of each neighborhood. 

We will be interested in the scaling limit in which we 
first let A ^ Z'^ (i.e. the systc ;m size tends to infinity) 
with X, y fixed, then let the Euclidean distance |x — y\ 
(in lattice units) become large. In this limit, the leading 
behavior of any of the ratios or probabilities just defined 
will behave as a power law, V oc \x — y\~'^'^'', where X/. is 
a scaling dimension (which also depends on d). However, 
in the antisymmetrized combination that appears in eq. 
(3), the leading part of 



J^ixivi,x2y2,-,xkyk) /J\f(x\x- y\ 



-2Xt 



(7) 



may be cancelled, so that a subleading power (with a 
larger exponent j^^"*'^y"™^ may be dominant. Thus in 
general the exponents must obey the inequalities 



j^antisymm ^ j^fc components ^ j^branches ^ j^paths /o^ 

As the addition of edges in the neighborhood of x 
or y contributes only a constant factor to the numbers 



M 



branches ; 



we have that X 



k components 
fe 



j^branches Yqx 



the two-dimensional case, the tree branches must enter 
the neighborhood of x in a sequence, and similarly at 
y. As the tree branches that cross from x to y cannot 
intersect because A is planar, the only nonzero terms in 
eq. (3) are for Ps that differ only by a cyclic permu- 
tation. For k odd, all cyclic permutations of k objects 
are even permutations, so all terms have the same sign, 
as if the sum were symmetrized. Therefore in d — 2, 

x^antisymm components r 7 jj ttt i j. j.i j. 

Aj, = Aj, ^ for k odd. We also note that 

in d = 2 dimensions, the boundary of a "thickened" tree 
is a nonintersecting dense loop, and the crossing by dis- 
tinct branches corresponds to crossing by the loop, 2k 
times. 

There is one further subtlety that must be mentioned. 
We have defined ratios and probabilities on finite graphs, 
followed by an infinite- volume limit. Each finite graph is 
spanned by a single tree, by construction, but as the vol- 
ume increases, a path from one vertex to another in a 
typical tree may involve larger and larger excursions, so 
that for infinite volume, two vertices might be connected 
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only "at infinity" , and then they can be regarded as not 
connected. Then the hmiting measure is for a forest of 
trees, not a single tree. It turns out that for d < A, 
there is a single tree in the infinite- volume limit, while 
for d > 4 there are infinitely many (infinite-size) trees [8] 
(these statements hold with probability one). Crossing 
probabilities by k distinct branches or k distinct paths 
can be defined here as well, with crossing not allowed 
to be "through infinity" ; we denote exponents under this 
condition by "< oo" . These probabilities will be less than 
or equal to the corresponding ones defined above, and 
so the related exponents obey x^"*"'^^^«<°° > Xj^^nches^ 

^paths<oo > ^paths j^^^ ^<oo = (rf _ 4)/2 for 

d> A (Ref. [8], Thm. 4.2), which because of the way this 
probability is normalized (any vertex is on some tree) is 
consistent with the belief that the Hausdorff dimension 
of any of the trees in the forest is 4 for d > 4. It would be 
interesting to obtain the remaining exponents x^'^^^^'^°° 
and ^^'■''"'='^°"<°° for A; > 1 also. In any case, the distinc- 
tion in definition disappears for d < 4 where there is a 
single spanning tree. 

We now turn to the calculation of the antisymmetrized 
combinations of Afs. According to eq. (3), the antisym- 
metrized sum of the relevant ratios is given by a ratio of 
cofactors: 

cof A^^i--^*'^!--^'^^ 
cof A(^i'?'i) 

PeSk 

As is well-known, such cofactors can be obtained from 
Gaussian integrals over Grassmann variables tp, tp*, 

cof A^^^'"^*"^^'"^*^^ 

where the overall sign is determined by the order of the 
Grassmann integrations (not by the selected values of Xi, 
Hi). In the limits A Z'^, followed by the limit of x 
and y far apart, the cofactors above become Gaussian 
integrals for a continuum massless complex scalar Fermi 
field. The equation of motion for the Fermi field "0 is 
simply A^" = 0, where A = Y^'^^^i 9ij.d^ is the Laplacian 
in d dimensions. As all Xi — > x, the ratio of cofactors 
becomes a sum of correlation functions of operators of 
the form 

o{x) = i,d^i>---d^,^,...i> (11) 

(with k tps, and where dfj,^^^ -- — '^Aii^M2 ' ' ') x, with 
a similar operator at y with ip* in place of tp. The 
undifferentiated ip and ip* are necessary to cancel the 
zero mode, both in the numerator and denominator of 
the correlation function. The remaining integrals can 
be simply expressed (using Wick's theorem) in terms 



of sums of products of derivatives of Green's functions 
G{x,y) = A~^[x,y) for the scalar field ijj in. d dimen- 
sions, and the required scaling limit of this expression 
exists; one has G{x,y) cx |x — y\'^~'^ for d > 2. Thus the 
scaling dimension oitp or ijj* is (d— 2)/2 in d dimensions, 
and an operator of the above form O has scaling dimen- 
sion ^^"t'^'y™™ = dimO equal to (fc — l){d — 2)/2 plus 
the number of partial derivatives in O (note that we re- 
placed A; by A; — 1 because the subtracted zero mode does 
not contribute to scaling). This implies that for fc = 1, 
the operator has dimension zero, which is correct as a 
spanning tree connects any two points x, y. It will be 
convenient to define dim'O = dimO — (fc — l){d— 2)/2. 

To find the operator that contributes the leading be- 
havior of the correlation function for a given fc, we must 
use as few derivatives as possible. Further, the equation 
of motion implies that any trace such as ^a'W3---V' 
vanishes. Then the leading term O is a product in which 
each multi-index partial derivative df_i^fi,^...ijj is a trace- 
less symmetric tensor, and the total degree (number of 
derivatives) in the product is as low as possible. Because 
of the anticommutation of the -08, O vanishes unless the 
traceless symmetric tensors d^.^^^,,,ip are linearly inde- 
pendent. We notice that the traceless symmetric tensors 
of given rank (degree) in dimension d, form an irreducible 
representation of the rotation group in d dimensions, 
SO(rf). In general, we expect that the leading part of the 
crossing probability Vk transforms as a scalar under ro- 
tations. A scalar (rotationally-invariant) operator can be 
obtained if the product includes either all or none of the 
members of a complete linearly-independent set of trace- 
less symmetric tensors for each degree (rank) /. Making 
the minimal choices of the degrees, we notice that this 
is analogous to filling states for fermions, where the sin- 
gle fermion states correspond to the symmetric traceless 
tensors. These tensors transform the same way as "hy- 
persphcrical harmonics" , which span the space of func- 
tions on a sphere S'^^^. This can be seen easily by rep- 
resenting each d/j, by a component of a vector x^, which 
transforms the same way, and the traces can be excluded 
if we assume that "^^Xfj^Xfj, is constant, so that tensors 
with non-vanishing trace are equivalent to lower-degree 
functions. Then the symmetric functions in .x^, binder 
this condition are simply the functions on S"^^^. If each 
fermion on 5"^"^ has a kinetic energy that is equal to the 
angular momentum /, then the many- fermion state with 
lowest total kinetic energy for given number of fermions 
fc corresponds to using the lowest total degree. The case 
where all traceless symmetric tensors of each (lowest) de- 
gree are used corresponds to filling a Fermi sea by filling 
the lowest shells up to angular momentum (degree of the 
traceless symmetric tensor) L. The total kinetic energy 
corresponds to the scaling dimension dim' O. When the 
lowest states are all filled, but the topmost shell is only 
partially filled, the scaling dimension interpolates linearly 
between the values it takes for filled shells. We point out 
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that the use of fermions on the sphere is more than an 
analogy, as the field theory of fermions on S"^^^ with time 
t corresponds to the original problem in radial quantiza- 
tion, obtained by conformalUy mapping R'^ to S"^^^ x R 
by a logarithmic change of variable, so that the dilatation 
operator becomes the Hamiltonian for the radial evolu- 
tion. 

We define N{1, d) to be the dimension of the space of 
traceless symmetric tensors of degree I in dimension d. 
When the shells are filled up to degree L, the preceding 
considerations lead immediately to the relations 



1=0 
L 

dim'Oz, =5I^^(^'^)- 



N{l,d) is given for all / > by 



N{l,d) = 



l + d-1 
I 



l + d~3 
I - 2 



(12) 
(13) 

(14) 



Here the first term is the number of symmetric tensors, 
and the subtraction is for removing the traces. Prom the 
binomial coefficients one sees that N(Ld.) is a polynomial 
in I of degree d — 2 for d > 2, and hence /c is a polynomial 
in L of degree d—1, and dim Ol is a polynomial of degree 
d. 

The leading behavior for / large is 



N{1, d) 



21 



d-2 



id -2)1 



(15) 



(throughout this paper, we use notation X ^ Y as Z ^ 
oo in the strict sense: lim2_>oo X/Y = 1). Then we find 



dim' Ol 



and hence 



dim'Oi ~ dimOi 



2L''-i 

d{d-2)\ ' 
-(d-1)! 



(16) 
(17) 

k^, (18) 



for the values of k specified. As mentioned above, for 
other values of fc, the scaling dimension (now for an oper- 
ator with nonzero spin in general) lies on a piecewise lin- 
ear continuous curve that interpolates the values above, 
and lies above that given implicitly by eqs. (12), (13) as 
polynomials in L, which are trivially extended to contin- 
uous values. 

By the inequalities (8), this result gives only an up- 



per bound on the leading exponents X\^ 



branches 



or X, 



paths 



However, the general formulas do give the exact expo- 
nents x™*'^^™™ for some, possibly subleading, terms in 



the probability. The rate of growth of the dimensions 
j^paths ^j^g ^^gg cross from x to y was shown rig- 
orously to be less than of order k'^/i'^~'^) in Ref. [9]. We 
obtain a bound with the same power, but now with a 
precise coefficient, and with subleading corrections. Note 
that the piecewise-linear curve for x^"*'»y™™ jg very close 
to its lower envelope, close enough that they have the 
same average rate of growth, eq. (18). 

We now consider the exact form of the dimensions ob- 
tained here for the k values given by eq. (12) in small d. 
In two dimensions, N{1, 2) — 2 {I > 0), and then 



dimOfc = (fc2 - l)/4 



(19) 



for k odd. This is in agreement with earlier results 
[2,3,5,6]. (Note that the scaling exponent for the path 

connecting x to y along the tree [4] corresponds to the 
case A; = 2, by considering the dual tree.) As emphasized 
above, for d = 2 and k odd, the arguments in this paper 
produce the exact ^ j^^branchcs ^ (^2 _ ^y^^ 

not only a bound. Note that after replacing k by k/2, 
this result is the same as the "fc-leg" dimension for k 
crossings by a dense polymer [10]. 

For d = 3, we have the familiar formula N{1, 3) = 21+1, 
and then k = {L + 1)'^, so L = \fk — 1. For the scaling 
dimensions, dim' Ol = L{L + 1)(4L + 5)/6, and 



-^antisymm ^ ^. ^ ^ 2^3/2 _ 1^1/2 
K 3 



(20) 



1 

6" 2' 

For d > 3, one can similarly solve explicitly, but the 

results are not as simple (in particular, they are not poly- 
nomials in k'^^). I am grateful to C. Tanguy for point- 
ing out that eqs. (12), (13) can be summed in closed 
form for all d, and that for d odd, fc is a polynomial in 
[L + {d— l)/2]^ of degree (d — l)/2. Hence in the cases 
d = 4, 5, 7, and 9, dimOfe can be expressed in terms of 
radicals in k. For d = 6, 8 and all d > 10, one meets 
polynomials of degree greater than four, and the results 
presumably cannot be expressed in terms of radicals. 

It is tempting to believe that the results obtained here 
for the "filled shell" values of k, and their smooth (poly- 
nomial in L) extension to general k, might be the ex- 
act values of x^^anches ^paths general dimension 
d > 2, as well as for d = 2. While it appears quite pos- 
sible that X^'-a^'^hes ^ ^paths general, it is not at all 

clear that they equal x'^'^^^^™^^^ . especially as the equal- 
ity that holds in two dimensions could be obtained from 
a simple argument that all cyclic permutations of an odd 
number of objects are even, an argument that definitely 
does not go through in d > 2. 

In conclusion, we have obtained, essentially rigorously, 
a precise upper bound on the exponents for k crossings of 
the uniform spanning tree on a finite graph in d dimen- 
sions, as well as some exact scaling dimensions in each 
dimension. However, we have not addressed the expo- 
nents in the uniform spanning forest which obtains on an 
infinite graph for d > 4. 
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